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Abstract 

Harmonic and Intermodulation distortions occur when a physical system is excited with a single 
or several frequencies and when the relationship between the input and output is non-linear. 
Working with non-linearities in the Frequency domain is not straightforward specially when the 
relationship between the input and output is not trivial. We outline the complete derivation of 
the Harmonic and Intermodulation distortions from basic principles to a general physical system. 
For illustration, the procedure is applied to the Single Mode laser diode where the relationship of 
input to output is non-trivial. The distortions terms are extracted directly from the Laser Diode 
rate equations and the method is tested by comparison to many results cited in the literature. 
This methodology is general enough to be applied to the extraction of distortion terms to any 
desired order in many physical systems in a general and systematic way. 
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I. INTRODUCTION 



Harmonic and Intermodulation distortions occur when a signal at a single frequency 
or a superposition of several signals with different frequencies propagate in a non-linear 
physical system. Historically it was first highlighted in the Radar and Radio Communi- 
cations industry because of the interest in understanding interference from other radars, 
jammers, other transmitters or modulators... Presently, it is pervading High Energy 
Physics particle detectors, wideband amplifiers, satellite communications |]T[ and other 
areas of science and technology. To cite some of the areas of interest, it occurs in many 
types of devices such as Mechanical Acoustical (Microphones for instance), Electronic 
and Microelectronic [|J, Microwave 0, Optical ||, Magnetic [0 and Superconducting ||. 

Parasitic frequency terms appear either at integer multiple of the base frequency 
(harmonic) or as a mixture of two or several multiples of base frequencies (intermodulation) 
when several base frequencies are used (as in modulation systems for instance). These 
terms can be either post-filtered or the signal can be pre-distorted in order to avoid the 
appearance of these unwanted terms. 

In this work, we introduce a general and systematic method to evaluate these terms from 
the equations of motion describing the physical system at hand (Laser diode rate equations). 
This case is chosen to highlight the case when the relationship between input and output is 
non-trivial. This paper is organised as follows. In section 2 we outline the general procedure 
for expressing the output from the input and define the various distortion terms. In section 
3 we apply the procedure in detail to the Single Mode Laser diode order by order up to third 
and conclude in section 4. The Appendix contains the full expressions of all terms derived 
in section 3. 

II. GENERAL PROCEDURE 

The Single Mode (SM) laser diode rate equations are written as: 
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M = r E N_^ + rg(N _ Nt){i _ eS)s (2) 

N represents the electron density ( N t at transparency) and S the photon density. r n 
is the electron spontaneous lifetime and r p h is the photon lifetime. (5 is the fraction of 
spontaneous emission coupled into the lasing mode, Y the optical confinement factor, g is 
the differential gain and e is the gain compression parameter, q is the electron charge, V 
the volume of the active region and I is the injection current. 

a- Step 1 [Elimination of (N - N t )]: 
The input is the injection current I and the output is the light intensity represented by 
S. Since other variable such as N, the electron density, intervene in the SM equations 
we proceed to the elimination of that variable. When several modes are present (or other 
intermediate variables) we proceed in an analogous manner by successive elimination until 
we have a single equation relating input to output. 



From (2) we extract the value of (N — N t ) as: 
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This can be used to find dN/dt and eliminate completely N from the coupled equations 
(1) and (2). The result is: 
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b- Step 2 [Small signal expansion about a static operating point]: 
The useful SM Laser diode regime with an injection current modulating the output light 
intensity. In the small dynamic signal case, I = Io + i and S — So + s, where Jo and So are 
the static injection current and photon density respectively, whereas, i = i(t) and s = s(t), 
are the signals of interest. We expand (4) [to n-th order]: 



i = {A lS + B lS ' + C lS "} + {A 2 s 2 + B 2 ss' + C 2 ss" + D 2 [sf}... + {A n s n 

+B n s n ~ l s l + C n s n - l s" + D n s n - 2 [s'f} (5) 

The expansion has this form because the first time derivative is applied to a fraction 
whose expansion contains n-order terms of the form s n ~ l s' and s n . After the derivation we 
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obtain n-order terms of the form s n_ V, s n ~ 2 [s'] 2 and s n ~V. s' and s" are first and second 
time derivatives of s. In addition we still have n-order terms of the generic s n and s's" 1 ' 1 
forms coming from the expansion of the right hand side of (4). Hence the general n-order 
term is A n s n + B n s n ~ l s' + C n s n ~ l s" + D n s n ~ 2 [s'} 2 where all the parameters A n , B n , C n , D n 
depend on the laser parameters. The appendix contains explicit expressions for some of 
these coefficients. The absence of a constant term means that no dynamic response % exists 
when there is no dynamic light excitation [s— 0]. 

Incidentally, we have used for simplicity of illustration the SM standard equations that 
contain two coupled population N (electron) and S (photon) equations in the form of two 
first-order ordinary differential equations. The method is in fact valid for any number of 
population equations as long as the elimination procedure of all intermediate variables (or 
poulations) is possible, leaving us with a single equation relating i to s (Equation (5)) or 
input to output. 

The harmonic distortions and intermodulation distortions are calculated from the generalized 
transfer functions denoted as H n (ui, 002,003...). They are obtained from the Fourier transform 
of the Volterra impulse response h n in the following way: 

/oo /*oo roo 
/ ... / h n (t 1 ,t 2 ,t 3 ,...)e- j ^ tM2+ ^ t3 ---Ut 1 dt 2 dt 3 .^ (6) 
-00 J — 00 J— 00 

For instance, the n-th order distortion is given by: 

J_- J_- ^ on 1 !\ H n(±JWl,±jW2,±jW3-)\\ / 7 n 

M„ (± „, ± «, ± «...) = 20 log,, 1 2 „- ln:Ui|gi(±J , m)l } m 

c- Step 3 [Method of Harmonic Input allowing the calculation of the Volterra transfer 
functions H n ]: 

The Harmonic Input Method (HIM) allows us to find directly all the if n 's in the following 

way: 

1. Express i as the sum : exp(ju>it) + exp(joo 2 t) + exp(jco 3 t)... 

2. Express s as the sum: 



00 

s = G k i m ..exp[j(kuJi + Iuj 2 + mw 3 ...)t] 



(8) 



The different H^s are found by direct identification of the Gkim... coefficients. For in- 
stance, we have: 

Gooo = 0,Gioo = #i(M),Giio = H 2 (ju 1 ,ju 2 ),... (9) 

III. ORDER BY ORDER DISTORTIONS 

Starting from the two laser rate equations, Darcie et al. derived formulae pertaining 
to the second order, third order and intermodulation distortions for a channel excited by 
a superposition of two signals with angular frequencies u)\ and u 2 ■ The second and third 
order distortions are calculated respectively at 2u\ and 3uj\ whereas the intermodulation 
distortions are evaluated at 2u\ — uj 2 and 2u 2 — uj\. From our formula (7), we can write 
these distortions as [0, HOI: 



a- Second order: 



M„(M,M) = 20 log 10 { !f; (i "'„ J 7 )l M I < 10 > 



b- Third order: 



w/. ■ \ on i [ #3 J^i,J^i,J^i ) nn 

M 3 {juujwuju! = 20 log 10 <^ , f 11 

c- Intermodulation: 

AjfC ■ • \ on 1 / l-^3(M, M, -jw 2 )| 1 - , 

M 3 (jwi, jwi, -j^ 2 = 20 log 10 ^ — — , f , 12 

■ ■ x on , [ 1-^3(^2,^2, ~M) I ) /iq\ 

M 3 [ju 2 ,ju 2 ,-ju l ) = 2Qlog 10 l ——. , \ 13 

The HIM allows us to calculate the values of the various H^s for n=l, 2 and 3. As an 

illustration of the procedure, we calculate Hi and H 2 after performing steps 1 and 2 and 
having obtained the expansion (5) to the specified order [§, [T0f| . First, we calculate ifi after 
simply using i = exp(jujt) in (5) and: 

00 

s = EW(M*l ( 14 ) 

k=0 
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The term G\ [with the identification of the terms multiplying exp(jut)] obeys the relation: 

1 = [A x + jujBx - cj 2 Ci]Gi (15) 

where A\,B\ and G\ depend on the laser parameters. Using the HIM [see (9)] we can 
write: 

Hi (juj) = G 1 = 1/[A 1 +jujB 1 - uj'd] (16) 



The modulus of H\{ju) equal to 1/y {A\ — u 2 Ci) 2 + (ooBi) 2 is what Darcie et al. || call 
the small-signal frequency response R(lu). In order to calculate second order terms, we use 
i(t) = exp(ju)it) + exp(juo 2 t) in equation (5) [truncated to second order] along with: 

oo 

s= Gkiexp\j{ku x + lu 2 )t] (17) 

k,l=0 



As stated in (9) the various Gki are obtained from the following || [I(J: G 00 = 0., G 



r io 



Hi(ju>i),Goi = Hi(ju 2 ) and Gn = Hi{ju)i,ju)2). Identification of the terms multiplying 
exp{j\uj\ + u)2\i) yields: 

= [Ai + j(u x + u 2 )B 1 - (wi + cj 2 ) 2 Ci]Gh + [A 2 G iG w + A 2 G 10 G 01 + juo^G^G^ 
+jcu 2 B 2 GioGoi — lu 2 C 2 GoiGiq — ulC 2 GioGoi — luiuj 2 D 2 GoiGio — ujiuj 2 D 2 GioGoi] 

(18) 

With the use of (16) this can be written as: 
H^jujt, juj 2 ) = Gn = 

-[2A 2 +j(u 1 + u 2 )B 2 - K 2 + u 2 2 )C 2 - 2uj l uj 2 D 2 ]H 1 {juj 1 )H 1 {juj 2 )H l {j[uj l + uj 2 }) (19) 

The second order distortions are obtained from (10) and (20) once the expressions of 
A 2 , B 2 ,C 2 , D 2 are found from the direct expansion of (4). The Taylor expansion to the 
third order is made in the Appendix allowing calculations of third order and intermodula- 
tion distortions. Also, the values of the various coefficients are given as functions of the laser 
parameters. 



IV. CONCLUSION 



In this work, we outline a general and systematic method to evaluate Harmonic and 
Intermodulation distortions that occur when a signal at a single frequency or a superposition 
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of several signals with different frequencies propagate in a non-linear physical system. The 
method is generalisable to any order and pertains to cases where a non-trivial relationship 
exists between input and output exists. 

In the case of a SM Laser diode, we have a system of two coupled differential population 
equations describing the system, nevertheless a general procedure has been described in 
order to relate the input injection current to the output light intensity. The results obtained 
in this paper agree with those published in the literature and the methodology highlighted 
provides a self-contained framework to evaluate distortions in a systematic way. 
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Procedure: 

In order to calculate the various A n , B n , C n and D n terms ( for n=l, 2, 3) a Taylor expansion 
about S = S and s' = is to be made to third order of equation (4). The following term 
should be expanded before taking its time derivative: 



ds _ rgjy, 

dt 



Tph 



(20) 



Tg(l-eS)S + ^ 

We can simplify the expansion considerably using the following orders of magnitude valid 
in most single-mode laser diodes: (5 = 0(~ 1CT 3 to 10~ 4 ), r n = oo , 1 — eSo is practically 1 
and gS l/r ph . The term becomes after these considerations: 



d,S | g 
dt r ph 



Tg(l-eS)S_ 
The next term to be expanded to third order is: 

ds _ rpm 



(21) 



1 



[g(l-eS)S + - 



dt 



Tg{l-eS)S + ^ 



(22) 



When we account for the aforementioned orders of magnitude it becomes: 



dS | S 
dt Tph 



(23) 



The third order expansion of (21) after taking the time derivative gives: 

[e/(T ph rg)}s> + [l/(r^ )K + [e/(rgS ) - l/(r^ 2 )]([ S f + ss") 
+ [e 2 /(r^ ) - (3eS - 2e 2 S 2 - l)/(r^ 3 )]( S V + 2s[sf) 



(24) 



Adding the contributions from (23) [s'/T + s/(T ph T)] and term by term identification 
gives for the various orders: 



First order: Ais + Bis' + C\s" where: 

A 1 = l/{r ph T) 

B 1 = e/(r ph Tg) + 1/r 

d = i/(r^ ) (25) 
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Second order: A 2 s 2 + B 2 ss' + C 2 ss" + D 2 [s'} 2 with: 

A 2 = 
5 2 = 

C 2 = [e/(TgS ) - l/(TgS 2 )} 

D 2 = C 2 (26) 

Third order: A 3 s 3 + B 3 s 2 s' + C 3 s 2 s" + D 3 s[sf with: 

A 3 = 
B 3 = 

C 3 = [e 2 /(TgS ) - (3eS - 2e 2 S 2 - l)/(TgS*)} 

D 3 = 2C 3 (27) 
The frequency response: 

\H x (Ju)\ = 1/ [A 1V /(1 - tfd/Atf + {uB 1 /A 1 )^ (28) 
and the second harmonic distortion given by (10) and transformed with the help of (19): 

M 2 (M,M) = 20 log 10 {2cu 1 2 |C 2 iJ 1 (2j^ 1 )|} (29) 

agree with the results of Darcie et al ||. 
In order to calculate third order and intermodulation effects, we use eq. (5) truncated to 
third order and the values of the A n , B n , C n coefficients (n=l, 2, 3): 

i = [Ais + B lS ' + C x s") + C 2 [ss" + [sf] + C 3 [s 2 s" + 2s[s'} 2 } (30) 

Apply the HIM to the above equation with: i = exp{joj\t) + exp(jcu 2 t) + exp(ju 3 t) and: 

oo 

s = G k i m exp\j(kuj 1 + lu 2 + mu 3 )t] (31) 

k,l,m=Q 

Proceeding as before (eq.18) we calculate the Gkim coefficients as well as the various 
Volterra transfer functions: 

Gooo = 0,Gioo = #iCM),Goio = Hi(ju 2 ), G QO i = Hi{ju 3 ), 
Gno = H 2 (j(v u juj 2 ),G m = H 2 (juj u juj 3 ) 

Con = H 2 (ju 2 , juj 3 ) and finally G 1U = H 3 (ju 1 , ju 2 , ju 3 ) (32) 



from eq.(30). 

Collecting the terms multiplying exp\j(u)\ + uj 2 + uj 3 )t) gives: 

= [At + j(u x + u 2 + uj 3 )B 1 - (ujx + u 2 + u; 3 ) 2 Ci]G'iii 

— 2C 2 {u\ + to 2 + ^3) 2 [2GiioGooi + 2G101G010 + 2GiooGmi] 

-3C 3 (^i + u 2 + cj 3 ) 2 [6G'iooGoioG'ooi] (33) 

Using the above relations between the Gkim and the H n 's and (16) we get: 

H^juj!, jw 2 , ju 3 ) = (uji +u 2 + u 3 ) 2 Hi(j[ui +u 2 + u 3 \)[AC 2 {H 2 {jui,ju) 2 )Hi{ju 3 ) 

+H 2 (ju 1 ,juj 3 )H 1 (juj 2 ) + H 2 (juj 2 , jw 3 )-^i(M)} 
+lSC 3 H 1 (ju; 1 )H 1 (ju; 2 )H 1 (ju 3 )} 

(34) 

Using expressions (16) for H\{juS) and (19) for H 2 (ju>i, ju 2 ) along with eq.(ll), (12) and 
(13) the third order and intermodulation distortions can be evaluated and they agree again 
with Darcie et al. || results. 



10 



